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Abstract. The superintegrability, wavefunctions and overlap coefficients of the 
Dunkl oscillator model in the plane were considered in the first part. Here finite- 
dimensional representations of the symmetry algebra of the system, called the 
Schwinger-Dunkl algebra sd(2), are investigated. The algebra sd(2) has six generators, 
including two involutions and a central element, and can be seen as a deformation of 
the Lie algebra u(2). Two of the symmetry generators, J3 and J2, are respectively 
associated to the separation of variables in Cartesian and polar coordinates. Using the 
parabosonic creation/annihilation operators, two bases for the representations of sd(2), 
the Cartesian and circular bases, are constructed. In the Cartesian basis, the operator 
J3 is diagonal and the operator Ji acts in a tridiagonal fashion. In the circular basis, 
the operator J2 is block upper-triangular with all blocks 2x2 and the operator J3 acts 
in a tridiagonal fashion. The expansion coefficients between the two bases are given by 
the Krawtchouk polynomials. In the general case, the eigenvectors of J2 in the circular 
basis are generated by the Heun polynomials and their components are expressed in 
terms of the para-Krawtchouk polynomials. In the fully isotropic case, the eigenvectors 
of J2 are generated by little —1 Jacobi or ordinary Jacobi polynomials. The basis in 
which the operator J2 is diagonal is then considered. In this basis, the defining relations 
of the Schwinger-Dunkl algebra imply that J3 acts in a block tridiagonal fashion with 
all blocks 2x2. The matrix elements of J3 in this basis are given explicitly. 
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1. Introduction 

This is the second part of this series concerned with the analysis of the isotropic Dunkl 
oscillator model. In part I, the model has been shown to be superintegrable, the 
wavefunctions have been obtained in Cartesian and polar coordinates and the overlap 
coefficients have been found [5j. In the present work, the representations of the symmetry 
algebra of the model, called the Schwinger-Dunkl algebra (see below), are investigated. 
As shall be seen, this study entails remarkable connections with special functions such 
as the Heun, little —1 Jacobi and para-Krawtchouk polynomials. 

1.1. Superintegrability 

One recalls that a quantum system defined by a Hamiltonian H in d dimensions 
is maximally superintegrable if it admits 2d — 1 algebraically independent symmetry 
generators Si that commute with the Hamiltonian 

[S u H} = 0, l<i^2d-l, 

where one of the symmetries is the Hamiltonian itself, e.g. Si = H. Moreover, a 
superintegrable system is said to be of order i if i is the maximal order of the symmetries 
Si in the momentum variables. 

1.2. The Dunkl oscillator model 

The isotropic Dunkl oscillator model [2, 5, llj in the plane is possibly the simplest two- 
dimensional system described by a Hamiltonian involving reflections. It is second-order 
superintegrable and is defined by the Hamiltonian [5] 

n = -\ [{V^f + [V^f] + 1 -[x 2 + y% (1) 

where T> x * i is the Dunkl derivative jU [H] 

V^=d Xi + ^(l-R Xi ), d Xi = ^-, 

with I denoting the identity operator and R Xi , Xi G {x,y}, standing for the reflection 
operator with respect to the plane Xi = 0. Hence the reflections R x , R y that appear in 
the Hamiltonian (CQ) have the action 

Rxf(x, y) = f(-x, y), Ryf(x, y) = f(x, -y), 

and thus evidently R x . = I. In connection with the nomenclature of the standard 
harmonic oscillator, the model is called isotropic because the quadratic potential is 
50(2) invariant. For the full Hamiltonian flTJ) to have this symmetry requires of course 
that jj, x = jji y . 
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The Schrodinger equation associated to T-L is separable in both Cartesian and polar 
coordinates. The spectrum of energies £ is given by 

£ N = N + /i x + Hy + l, N = n x + n y , (2) 

where n x , n y are non-negative integers. The wavef unctions are well defined for the values 
H X} fly G (— |, oo); the = corresponds to the standard quantum harmonic 

oscillator. It is easily seen from fl2]) that the energy level En exhibits a iV + 1-fold 
degeneracy. 



1.3. Symmetries of the Dunkl oscillator 

The symmetries of the Dunkl oscillator Hamiltonian ([T]) can be obtained by the 
Schwinger construction using the parabosonic creation/annihilation operators [5| |8| ITU]. 
We consider the operators [T^l 115] 

A£ = -j=( Xi ±V£?), x^^y}. (3) 

It is verified that the operators A± satisfy the following commutation relations: 

[A X J, A°£] = I + 2fi Xi R xv {A±, R Xt } = 0, (4) 

where {x, y} = xy + yx denotes the anticommutator. In addition to the commutation 
relations ()4j), one has 

[A£ , Jg] = [A^ , R Xj ] = [R Xi ,R Xj ] = 0, i± j. (5) 

In terms of the operators ([2]), the Hamiltonian ([T]) takes the form 

U = \{A% A*_) + 1 -{A\, Al} = n x + H y , 

where 

n xi = \{A^,A x J} = -\{V^f + l -xl (6) 

is the Hamiltonian of the one- dimensional Dunkl oscillator. 

The symmetry generators of the Dunkl oscillator model are as follows. Consider 
the operator 

J* = \{A%A%}-\{At,A\}. (7) 
It is directly verified that [H, J3] = 0. Since J3 can be written as 



J3 — ~ 
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using (E]), it is clear that this symmetry corresponds to the separability of the 
Schrodinger equation in Cartesian coordinates j5]. A second symmetry generator is 
given by 

j 2 = ^{A%A y _-A°LA\). (8) 

It is again directly verified that [ J 2 , H\ = 0. In terms of Dunkl derivatives, this operator 
has the expression 

■h = \ % (xV^y - yV»*) , 

and it has been shown |5] that J 2 is the symmetry corresponding to the separation of 
variables in polar coordinates. A third symmetry J\, algebraically dependent of J 2 , J3, 
is obtained by taking J\ = — i[J 2 , J3]. This additional symmetry generator reads 

J x = l(AlA y _ + A°LAl). (9) 

In addition to Jj, i = 1, ... ,3, it is directly checked that the reflections R x , R y also 
commute with H. 



1.4- The main object: the Schwinger-Dunkl algebra sd(2) 

The symmetry algebra of the Dunkl oscillator, called the Schwinger-Dunkl algebra, is 
denoted sd(2) and defined by the commutation relations 

{J u R Xi } = 0, {J 2 ,R x .} = 0, [J 3 ,R Xi ) = 0, (10a) 
[J2,Js]=iJi, [J3,Ji]=iJ2, (10b) 

[Jl, h] =i\Jz + J^xRx + l^yRy) ~ U(HxRx ~ fJLyRy)/2\, (10c) 

where R x . = I, Xi G {x, y}, and the Hamiltonian % is a central element. The algebra 
sd(2) admits the Casimir operator [7] 

C = Ji + J 2 + J 3 + ij, x Rx/2 + /jyRy/2 + fix/^yRxRy, 

which commutes with all the generators. In the present realization, the Casimir operator 
C takes the value 

C= 1 -(H 2 -1). 

Note that the involution P = R x R y also commutes with all the generators and thus 
can be viewed as a second Casimir operator. Furthermore, it is easily seen that when 
— l^y — 0, the Schwinger-Dunkl algebra sd(2) reduces to the Lie algebra u(2), which 
is the symmetry algebra of the standard isotropic 2D quantum oscillator in the plane. 

The irreducible representations of sd(2) can be used to account for the degeneracies 
in the spectrum of H. In finite-dimensional representations of degree N + l, the action of 
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the symmetry generators J 1; J 2 , J3, -R^ and R y indicate how the degenerate eigenstates 
of % corresponding to the energy value Sn transform into one another under the action 
of the symmetries. In the following, three bases for the finite-dimensional irreducible 
representations of sd(2) will be constructed and the explicit formulas for the action of 
the symmetry generators on each basis will be derived. 

1.5. Outline 

Here is the outline of the paper. In Section 2, we construct the Cartesian basis in 
which the symmetry generator J3 is diagonal and J 2 acts in a tridiagonal fashion. 
In Section 3, we introduce the circular creation/annihilation operators and study the 
associated circular basis in which J 2 is block upper-triangular and J3 is tridiagonal. We 
show that the interbasis expansion coefficients involve the Krawtchouk polynomials and 
we derive the spectrum of J 2 algebraically. In Section 4, we obtain the eigenvectors 
of J 2 in the circular basis for odd-dimensional representations and show that these 
eigenvectors are generated by the Heun polynomials and that their components are 
para-Krawtchouk polynomials. The fully isotropic case is shown to involve the little —1 
Jacobi polynomials. In Section 5, the eigenvectors of J 2 in the circular basis for even- 
dimensional representations are studied. In Section 6, we examine the basis in which J 2 
is diagonal and show that J 3 acts in a six-diagonal fashion on this basis. We conclude 
with an outlook. 

2. The Cartesian basis 

In this section the Cartesian basis for the finite-dimensional representations of sd(2) is 
constructed using the realization (EJ), (ED, (EJ) of the algebra generators in terms of the 
creation/annihilation operators fl3]). The representation spaces spanned by the Cartesian 
basis correspond to the spaces of degenerate wavefunctions with energies £n,N £ N, 
separated in Cartesian coordinates, although a different normalization is used for the 
basis vectors. The action of the sd(2) generators on the wavefunctions were obtained by 
a direct computation in [5] using the expressions of the symmetries in terms of Dunkl 
derivatives. Here the actions of the generators and the spectra of the Hamiltonian 7i 
and the symmetry generator J3 are obtained in a purely algebraic manner. 

The Cartesian basis vectors are labeled by two non- negative integers n x , n y and are 
denoted by \n x ,n y ). These basis vectors are defined by 

\n x ,n y ) = (A x + ) n °(Al) n y\Q x ,Q y ), (11) 

where 10^,0^) is the "vacuum" vector. The vacuum vector has the defining properties 

A*|0 X ,0„)=0, A y _\0 x ,0 y ) =0, (12a) 

R X \0 X , 0y) = \0 X , 0y), Ry\0 X , 0y) = \0 X , y) . (l2b) 
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The action of the reflection operators and the creation/annihilation operators on 
the Cartesian basis vectors can be derived from the above definitions and the 
commutation/anticommutation relations (j3J) and (JSJ. From the anticommutation 
relations 

{Al,R x } = 0, {A v + ,R y } = 0, 

and the vacuum parity conditions (j!2b|) . it directly follows that 

R x \n x ,n y ) = (-l) nx \n x ,n y ), R y \n x , n y ) = (-l) ny \n x , n y ). (13) 

By the definition of the basis vectors (ITT]) , one has also 

A x + \n x ,n y ) = \n x + l,n y ), A y + \n x , n y ) = \n x , n y + 1). (14) 

To derive the action of the operators A x _i on the Cartesian basis, one needs the 
commutator identity 

[A x ~, (4?)"] = (A^ 1 [n + ^,(1 - (-l) n )R Xi ] , (15) 
which is easily proven by induction. It is convenient to introduce the //-numbers [H] 

[n]„ = n + Ml -(-!)")• (16) 
Using the identity (1T5]) and the formulas (ITTjl . (I12aj) and (112bl) . one finds 

A x _\n x ,n y ) = {A\) n y[A x _,{A%) n *]\U x ,%) = [n^Jn, - 1, n y ), (17) 

and similarly for A v _. 

Using the results (TLB"]) . (TT4"|) and (fT7|) along with the expressions of the symmetry 
generators Jj, i — 1, 2, 3, and in terms of the operators A± given in (JTJ), (JH]) and ([9]), 
one finds that the action of the symmetries on the Cartesian basis is given by 

J2\n x ,n y ) = ^([nJ^K + 1,% - 1} - [n x ]^\n x -l,n y + l)j, (18a) 

Mn x ,n y ) = - (n x - n y + /i x - fi y ) \n x ,n y ), (18b) 

and the action of J\ can be obtained directly by commuting J2 and J3. The central 
element % has the action 

U\n x ,n y ) = (n x + n y + \i x + // y + 1)1^,^). 

Hence the spectra of the symmetry generator J 3 and of the full Hamiltonian % of the 
Dunkl oscillator have been recovered in a purely algebraic manner. As is expected, the 
symmetry operators Ji, . . . J3 and the involutions R x , R y transform the set of vectors 
\n x ,n y ) with a given value of iV = n x + n y into one another; these vectors are the 
degenerate eigenvectors of % with energy Sn = N + fi x + \i y + 1. 
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The preceding results can be used to define an infinite family of N + 1-dimensional 
irreducible modules of the Schwinger-Dunkl algebra sd(2) (fTUj) . Let /j, x , fi y G R be real 
numbers such that fi x , fi y G (—1/2, oo) and denote by vjf x the iV + 1-dimensional C- 
vector space spanned by the basis vectors Vn , nS {0, . . . , N}. Consider the vector 
space vj^ x endowed with the following actions of the sd{2) generators: 

Jiv^ = \([N- nl y v^ + [nUv { :-f v) ) , (19a) 

J2v { r^ = Yi( [N ~ n ^ v ^ Ly) - N^tr y) ) , (wb) 

JA m) = U+ \{fi x -fiy- N)j (19c) 

where [n] denotes the /z-numbers ffT6j) . The central element % and the Casimir operator 
have the actions 



and 



Cv^ = \{(N + ii, + to){N + 2 + ii x + liy)}v^\ 

It is clear that yj^ 1 *'^ i s a srf(2)-module and its irreducibility follows from the fact 
that the /i- numbers appearing in the matrix elements of J 1; J 2 are never zero for 
fi x ,fi y G (— 1/2, oo). For /i x = fi y = 0, it is directly seen that the s<i(2)-module V^ ' ^ 
reduces to the standard iV + 1-dimensional irreducible su(2) module. 

3. The circular basis 

In this section, the circular basis for the finite-dimensional representations of sd(2) is 
constructed using the left/right circular operators. The actions of the symmetries on this 
basis are obtained and the spectrum of the generator J 2 is derived from these actions. 
The expansion coefficients between the circular and Cartesian bases, which involve the 
Krawtchouk polynomials, are also examined. 

The left/right circular operators for the 2D Dunkl oscillator are introduced following 
the analogous construction in the standard 2D harmonic oscillator [TJ. We define 

A L ± = -L (Al =f iAl) , A* = -±= {A% ± iAl) , (20) 

where A± are the creation/annihilation operators of the Dunkl oscillator that obey the 
commutation relations (j3J). The inverse relations are easily seen to be 

Al = ^=(A L ± + Al), Al = ^L(A L ± -Al). 
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The left /right operators obey the commutation relations 

[^,^3 = 0, [A$.,A*\ = 0, 

y ' 



[A^, A + ] — /Ji X R x ~ l^yRyj [ A -1 A^] — fl x R x — HyRy 



[A L _, A%\ = I + ii x R x + HyRy, [A*, A*} = I + fi x R x + fi y R y , 

and the algebraic relations involving the reflections become 

R X A± = —A±R X , R X A± = —A±R X , R y A± = A±R y , R y A± = A±R y . (21) 

The circular basis vectors \riL,n R ) are labeled by the two non-negative integers rii, ur 
and are defined by 

\n L ,n R ) = (A^(A*) n *\0 L ,0 R ), (22) 

where |0l,0.r) is the circular vacuum vector with the properties 

A L \0 L ,0 B ) = 0, A*\0 L ,0 R ) = 0, (23a) 
R X \0 L , Or) = \0 L , R ), R y \0 L , Or) = \0 L , Or). (23b) 

Given the definition (122|) . one has 

A+\n L ,n R ) = \n L + l,n R ), A+\n L ,n R ) = \n L ,n R + 1). 

From the relations (12T|) and the definition (1221) . it follows that 

Rx\n L ,n R ) = (-l) nL+n *\n R ,n L ), R y \n L} n R ) = \n R ,n L ). (24) 

Consider the commutator identities 

n 

[At (A L + r +i ] = (n+ i)(A L + r + j2( AL + r~ a «r {(-i)>*^ + /w , 

n 

[At (A*)^ 1 ] = J2( AL +) n ~ P Kf {(-1)"-^ - VyRy} , 

(3=0 

which can be proven straightforwardly by induction. From the definition (122)) . the 
vacuum conditions ()23ap . (123bp and the above identities, the action of A 1 ^ on the circular 
basis elements \ni,n r ) can be derived by a direct computation. For til = n R , one has 

Atn L ,n R ) = n L \n L - l,n R ), 

For ul > ur, one finds 

TlL — l 

Atn L ,n R ) = n L \n L - l,n R ) + ^2 {(-l) nR+: > fi x + ^[ul + n R - j - 
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Finally, for ni < n R , one obtains 

A L _\n L ,n R ) = n L \n L - l,n R ) - {(-l)"^- 7 /^ + fjLy}\n L + n R - j - 1, j). 

3=n L 

To obtain the corresponding formulas for the action of A R , one needs the identities 

n 

[At(A*r +i ] = (n+i)(A^r + Y,(^r~ a (A^r{(-ir-y x R x +fi y R y }, 

n 

[At (A L + ) n+1 ] = ^i)^^ {(-l)V^ - HyRy) . 

/3=0 

Using the same procedure as for A L , we obtain the action of At. For til = n R , we have 

A R \n L ,n R ) = n R \n L ,n R - 1). 

When ni > n R , one finds 

n L -l 

A R \n L ,n R ) = n R \n L ,n R - 1) + ^ {(-l) nR+J /i x - /i y }\n L + n R - j - 

j=n R 

and for n R > Ul, the result is 

n R -l 

A R \n L ,n R ) = n R \n L ,n R - 1) - ^ {(-l) nfl+ > x . - /^}|n L + n fi - j - 1, j). 

F I Ft 

As is seen from the formulas, the operators A_ have the effect of sending the circular 
basis vectors \nL,n R ) to all circular basis vectors \il3r) with ii + j R = ni + n R — 1. 

In terms of the circular operators ( l20l . the symmetry generators and the central 
element % take the rather symmetric form 

Ji = l ({A L + , A R } - {At AD) , J 2 = \ ({At A R } - {A L _, A%}) , (25a) 
J 3 = \ ({AtA R } + {A^, A*}) , ft = i {{At A L + } + {At A R }) . (25b) 

F I Ft 

Using the above formulas and the actions of the circular operators A± , the matrix 
elements of the sd(2) generators in the circular basis can be computed; they are given 
below for J 2 and J3. The action of the Hamiltonian H is 

^L\n L ,n R ) = {n L + n R + fj, x + fj, y + l)\n L ,n R ). 

It is clear that the generators preserve the subspace spanned by the basis vectors 
{\ n Li n R) I n L + n R — N}. As is seen from the action of H, this corresponds to the 
space of degenerate eigenstates of H, with energy The properties of representations 
of the symmetry generators in the circular basis will now be used to derive the transition 
matrix from the circular basis to the Cartesian basis and to obtain the eigenvalues of 
J2 algebraically. 
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3.1. Transition matrix from the circular to the Cartesian basis 

We consider the N + 1-dimerisional energy eigenspace spanned by the circular basis 
vectors \til, Ur) with til + tir = N and redefine the basis vectors as follows 

Bi := {/„ = \0 L , N r ), h = \U, (N - 1) R ), ...,f N = \N L , B )}. 
On this basis, a direct computation shows that the generator J 3 has the action 

hU = 1{(N- n)f n+1 + Un + n/ n _i}, (26) 

where we have defined 

Since J 3 is diagonal in the Cartesian basis and tridiagonal in the circular basis, the two 
bases are related by a similarity transformation involving orthogonal polynomials. 

Let us consider the decomposition of the Cartesian basis vector y[^ x '^ of vjf x 
on the circular basis 

TV 

=J2 C nU)fn, (27) 
n=0 

where j G {0, . . . , N}. Acting on both sides of (I27p with J 3 and using fU9cp and ( I26p . 
one arrives at the following recurrence relation satisfied by the expansion coefficients 
C n (j): 

(2j - N)C n (j) = (n + l)C n+1 (j) + (N-n + l)C n ^{j), 
with C_i = 0. Upon factoring out the initial value 

Cn(j) = C (j)P n (j), 

we obtain the recurrence relation 

(2j - N)P n (j) = (n + 1) P n+1 (j) + (N — n + 1) P n -i(j), (28) 

where Po(j) = 1. It follows from (128 p that P n (x) is a polynomial of degree n in x. Upon 
substituting P n (j) = P n (j)/n\, we obtain the normalized recurrence relation 

(j - N/2)P n (j) = P n+l (j) + -n(N -n+ l)P„-i(i). 

It is directly seen that the polynomials P n (j) are the monic Krawtchouk polynomials 
K n (x;p,N) [H] with parameter p = 1/2 and variable x evaluated at x = j. We thus 
have 

C n (j) = C (j)K n (j;l/2,N), 

where the constant C*o(j) can be chosen to ensure the unitarity of the transition 
matrix by using the orthogonality relation of the Krawtchouk polynomials. Despite 
the differences that the Dunkl and standard harmonic oscillators exhibit, the relations 
between the circular and Cartesian bases are identical in both cases. 
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3.2. Matrix representation of J2 and spectrum 

The circular representation space can be used to derive the spectrum of the symmetry 
operator J 2 . To exhibit the structure of J 2 , we introduce the following notation for the 
basis vectors: 



\n L ,n R ) = \£,±), 

where 



I = U n L ~n R \/2\, ± = sign (n R -n L ), 
where \_x\ is the floor function. We adopt the convention that 

signO = — 1 

for convenience. We denote by B2 the circular basis such that + n R = N with the 
ordering 

£ 2 = {|0,+),|0,-),|1,+),|1, -),...,} 
As an example, consider the case N = 3. The basis reads 

£ 2 = {|0,+),|0,->,|l,+),|l,->} 
and corresponds to the following ordering of the standard circular basis vectors \ni, n R ): 

£ 2 = {|1,2>,|2,1>,|0,3>,|3,0>}. 

For iV = 4, one has 

i3 2 = {|0,-),|l,+),|l,-),|2,+),|2,-)}, 

which corresponds to 

B 2 = {|2,2),|1,3),|3,1),|0,4),|4,0)}. 



Using the action of the operators A± R and the formulas (125 aft . ( I25b[) . the matrix 
representation of J 2 in the circular basis i3 2 is derived in a straightforward manner. 
We find that for N even, the iV + 1-dimensional square matrix representing J 2 in the 
basis B2 is block upper-triangular with all blocks 2 x 2 in addition to a row of 1 x 2 
blocks. The matrix reads 



[MB2 



/0 (jJi Ll>2 ••• 

ri ni n 2 
r 2 Ox 



Mm \ 
^•m— 1 

T m _i ill 

-Trfi / 



(29) 
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with m = N/2 and where we have 



k + C/2 -C/2 
C/2 -fc-C/2 





(-a 





< 








f; 











odd, 



even, 



(30) 



with Uk corresponding to the lower part of Qfc. We have taken 



C = fix + fiy 

In the iV odd case, one obtains 

/To 



2j% 



fl 



? 2 



n 



r 



m— 1 



ra—1 



(31) 



with m = (JV — l)/2 and where 



(2fe + 1 + 0/2 
-C/2 



C/2 

-(2fc + 1 + 0/2 



A: odd, 



even, 



(32) 



Since in both cases the matrices representing J2 are block upper-triangular, it follows 
from elementary linear algebra that the set of eigenvalues of J2 is the union of the sets 
of eigenvalues of each diagonal block or IV By the direct diagonalization of the 2x2 
diagonal blocks, we obtain that when iV is even, the eigenvalues of J 2 are given by 



A 



± 



±Jk(k + fi x + fi y ), k = 0, 



m. 



where m = N/2 and where the eigenvalue A = is non-degenerate. When iV is odd, 
the spectrum of J 2 has the form 



± 



y / \k + fi x + l/2)(fc + fiy + 1/2), 



0. 



m 



where w! = (N — l)/2. These eigenvalues are indeed the eigenvalues of J 2 that 
were obtained in the first part |5J by solving the differential equation arising from the 
realization of J 2 in terms of Dunkl derivatives; here they have been obtained in a purely 
algebraic manner. It is seen from the matrices (j29p . (1ST]) that when fi x = fi y = 0, the 
matrices representing J 2 are diagonal. This corresponds to the standard result for the 
harmonic oscillator, where the circular basis is the eigenbasis of the symmetry J 2 . 

Given that in the case of the Schwinger-Dunkl algebra sd(2), the circular basis does 
not diagonalize J 2 directly, it is of interest to inquire about the eigenvectors of J 2 in this 
basis. This is the subject of the next two sections. 
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4. Diagonalization of J 2 : the iV even case 

This section is devoted to the computation of the eigenvectors of J2 in the circular basis 
for the N even case. To perform the calculation, we shall make use of an auxiliary 
operator Q whose eigenvectors have been related to those of J2 in the previous paper 
[5j. The evaluation of the eigenvectors of Q is somewhat involved and consequently it 
is instructive to first expose the main steps of the computation. 

Firstly, the structure of Q will be used to reduce the eigenvalue problem to 
a system of recurrence relations for the components of the eigenvectors. Secondly, 
using generating functions, the recurrence system will be transformed into a system of 
differential equations and the solutions will be expressed in terms of Heun polynomials. 
Thirdly, using well-known properties of Heun functions, the explicit expressions for 
the components of the eigenvectors will be obtained in terms of a special case 
of complementary Bannai-Ito polynomials [6J which correspond to para-Krawtchouk 
polynomials |18| . Lastly, the relation between the eigenvectors of J 2 and Q obtained in 
the first paper |5j will be used to write the final expression for the eigenvectors of J 2 in 
the circular basis. 



4-1. The operator Q and its simultaneous eigenvalue equation 



The operator Q has been used in the first part [5] to obtain the overlap coefficients 
between the wavefunctions in Cartesian and polar coordinates of the 2D Dunkl oscillator 
|7] . It is defined in terms of J 2 through the relation 



Q = —2iJ 2 R x — H-xRy — H>yRx — (l/2)R x R y . 



(33) 



Given the action (124]) of the reflections operators, it is seen that R x , R y have the following 
matrix representation in the circular basis £> 2 : 



R x = R y = diag( l,o-i,cxi, ■ ■ ■ ,a 1 ), a 1 




(34) 



Using the formula (I34p for the reflections and the formulas ( |29|) . (130]) for the expression 
of J 2 in the N even case, one obtains from ( 133]) 



/ ¥0 h S 2 
$1 Ai 



V 



m—1 



A m _l 



Ai 



with m = N/2, (p = — ( — 1/2 and where 



»C - 1/2 
2ik - (1 - i)C 



-2ik - + 
-iC-1/2 



t 


(-2ii 


2iA 




\-2ii 


2ii) 




(2K 


-2i(\ 




\2K 


-2i() 



modd 



meven 



(35) 
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The 1x2 blocks 5i correspond to the lower part of the blocks A,. From the block 
upper-triangular structure, it follows that the eigenvalues of Q are given by 

2fc + C-l/2, Pk = -(2& + C + V2), k = l,...,m, (36) 



and we also have v§ = —( — 1/2. Let us denote by \k, ±)q the eigenvectors of Q with 
eigenvalues v^. We wish to evaluate the components of these eigenvectors in the circular 
basis. We define their expansion in the circular basis by 



(37) 



£=0 
<r=± 



^=0 

cr=± 



for k = 1, . . . , m and where the vectors \£, ±) are vectors of the circular basis B^- It is 
clear from the matrix representation of Q that |0, — }q = |0, — ) and thus f^"(0) = 1. 

We shall study the simultaneous eigenvalue equation for the operator Q. Since the 
matrix representing Q is block upper-triangular, the matrix of eigenvectors will have 
the same structure. We define the matrix of eigenvectors of Q as follows: 



/I Vbi V 02 



W 



Vim 



22 



Kn-lm 
VrTn.ni. 



where 



Vt 



£,k 



u+(k) v+{k) 
uj{k) vj{k) 



and where Vtk is the 1x2 block corresponding to the lower part of V^. The simultaneous 
eigenvalue equation for the matrix [Q]b 2 can be written as 

W ■ L = [Q] B2 ■ W, (38) 

with 



L = diag(i/^ } , Ai,--- ,A m ), A k 



2fc + C - 1/2 

-2fc-C-l/2, 

As will be seen, the components (137|) of the eigenvectors of Q can be derived from the 
eigenvalue equation ( 138]) by solving the associated system of recurrence relations. 



4-2. Recurrence relations 

It will prove convenient to consider the two sectors corresponding to the eigenvalues v 



k ■ 



v k separately. In block form, for 



, m and k 



(138]) can be written in the form 

V lk k k = $ e V ek + E A i V i k > 



k-l 



, m, the eigenvalue equation 



(39) 



with $^ and Aj given in (135]) and where the range of the sum has been determined by 
the structure of the eigenvector matrix W. 
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4-2.1. The u£ eigenvalue sector 

We consider the eigenvectors \k, +}q of Q with the expansion 

k 

\k,+) Q = J2 u e(k)%°), (40) 
e=o 

<T=± 

and associated to the eigenvalue u£ — 2k + £ — 1/2. It is understood that Uq (k) does 
not belong to this decomposition. For £ — 1, . . . , m, it directly seen that the eigenvalue 
equation fl39|) is equivalent to the following system of recurrence relations: 

k 

[2k + (1 - i)(]u+ = [-2U - (1 + i)t\uj -2iJ2 {(-l)^Vx + Hy}Bj, (41a) 

j=e+i 

k 

[2k + (1 + i)£]uj = [2i£ - (1 - i)C]u+ - 2i {(-l) J "Vx + Vy}Bj, (41b) 



where we have defined 



B j = Uj - u+ 



and where the explicit dependence of the components uf on k has been dropped for 
notational convenience. The case £ = is treated below. The system of recurrence 
relations (|4ip is "reversed": the values of are obtained from the values of with 
i < j and j < p. The terminating conditions are at £ = k. In this case ( |4ll reduces to 



[2A: + (1 - z)C]< = [-(2A:)i - (1 + i)C]% , (42a) 
[2k + (1 + z)C]w^ = [(2k)i - (1 - (42b) 

In accordance to the system ( 1421) . we choose the following terminating conditions 

M fe = u k = L 

Upon introducing 

= w+ + uj, 

the system (T4T|) is directly seen to be equivalent to 

k 

[k + (]A e = -i{£ + C)B t -2iJ2 + (43a) 

j=£+l 

[k]B e = i£A e . (43b) 

The above system accounts for the £ = case. Indeed, it is seen that Bq = and 
hence uj = Aq/2. These equations can be simplified by factoring out the terminating 
conditions 

A e = a A e , B e = f3 B e , 
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where a = (1 — i) and /3 = (1 + i). It is seen that the normalized components Ag, B% 
are real and satisfy the system 

k 

(k + ()A e =[£ + (}B £ + 2 {(-^'"Va. + ( 44a ) 

j=M-i 

fci^ = M/, (44b) 

with the terminating conditions = B^ = 1- The system (144")) can be simplified by 
introducing the reversed components ae — ^fc-£ and 6^ = B^-^- Using the index n, the 
system takes the usual form 

n-l 

(k + CK = (k - n + OK + 2 ^{(-1)"+^ + N }b a (45a) 

kb n = (k — n)a n , (45b) 

with the initial conditions a = 6 = 1. Hence the components uf(k) of the eigenvector 
\k, ±)q of the operator Q have the expression 

U e (K) — - , U £ (KJ — - , (40 J 



where a n and 6 n are the unique solutions to the system (1451) . 
4-2.2. The eigenvalue sector 

We consider the eigenvectors \k, — )q corresponding to the eigenvalue of Q with the 
circular basis expansion 



£=0 
er=± 



and associated eigenvalue = —2k — ( — 1/2. An analysis similar to the preceding 
one shows that the components vf{k) differ from the components uf(k) only by their 
terminating conditions. Again choosing v^(k) = 1, we find 



This yields 



where 



lo^k-i + toh-e loa-k-e + £oh-i 



vj{k)= — ' — , ,+ (fc)= — ' — , (47) 



2i(fc + C) 2A; 



[l + i)k + iC (l + i)k + iC 

and where t>^ = ■joOk- 
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4-3. Generating function and Heun polynomials 

We have seen that the evaluation of the components of the eigenvectors of Q in the 
circular basis depends on the solution of the recurrence system ( 1451) . As it turns out, 
an explicit solution for a n (k) and b n (k) can be obtained using generating functions. 
We introduce the ordinary generating functions 

We shall make use of the elementary identities 

zd z A(z) = na n z n , (1 - z)' 1 ^) = £ I £ a k J z n , (48a) 

(1 + z)- l A(z) = £ I E (" 1 ) fc+n ^ I * n - ( 48b ) 



Using the above identities, it is easily seen that the system of recurrence relations (|45j) 
for the quantities a n , b n is equivalent to the following system of differential equations 
for the generating functions A{z), B(z): 

(k + Q)A(z) = (k-C- zd z )B(z) + Y^~ z B{z) + Y^Hz), (49a) 

kB{z) = {k- zd z )A(z). (49b) 

By direct substitution, we find that the generating function A(z) satisfies the second- 
order differential equation 

#(,) + (iz*Z< + *<„ + ( T^.+ fJ = 0. ,50) 

V z z-l l + zj \z(z - l)(z + I) J 

This corresponds to Heun's differential equation [3J d3j . The general form of the Heun 
differential equation is 

w"{z) + + —) w'(z) + aPz - q w(z) = 0, (51) 

\z z — 1 z — aj z{z — l)(z — a) 

with a + (3 + l = r y + 5 + e. Comparing f lBTj) with (15D|) . we thus write 

A(z) = H£(a,q;a,(3,j,5,z) (52) 

with the parameters 

a = — 1, q = 2k(fi y — fi x ), a = —2k, (53a) 
/3 = fi x + fly, 7 = 1 - 2k - fi x - fly, 5 = 2fi y . (53b) 

The function H£(a,q;a, (3,^,8) denotes the solution to (I5ip that corresponds to the 
exponent at z = and assumes the value 1 at that point. This is obviously the case of 
A(z). It will be seen that A(z) is in fact a polynomial of degree 2k, and hence that the 
Heun function (152]) is in fact a Heun polynomial. Given the system (|49|) . we also have 

B(z) = k~\k- zd z )A(z). 
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4-4- Expansion of Heun polynomials in the complementary Bannai-Ito polynomials 

The well-studied properties of Heun functions can be used to obtain a closed form 
formula for the coefficients a n and hence for b n . In what follows, it will be shown that the 
Heun polynomial in (152]) can be expanded in terms of a special case of the complementary 
Bannai-Ito polynomials corresponding to the para-Krawtchouk polynomials. 

Consider the solution H£(a,q; a,/?, 7, 5) to the equation fl5T]) and its Maclaurin 
expansion 

oo 

H£(a, q; a, (3, 7, S) = ^ c nZ n , 

n=0 

where c_i = 0, Cq = 1. The coefficients c n obey the three-term recurrence relation [3|1T3] 

Rn c n+l — {Qn + ?) c n + Pn c n-l = 0) 

where 

R n = a(n + l)(n + 7), Q n = n[{n - 1 + + a) + a5 + e] , (54a) 
P n = (n-l + a)(n-l + P). (54b) 

The identification of A(z) as a Heun function enables one to reduce the evaluation of 
a n to the solution of a three-term recurrence relation. It is seen that with the choice 
of paramaters (153]) . the expansion coefficients c n of the Heun function A(z) truncate at 
degree k = 2n. Hence A(z) is a polynomial of degree 2k in z. For convenience, we use 
the symbol V n = a n in the computations to follow. Using the parameters f l53|) in the 
recurrence coefficients (1541) for the expansion coefficients in 

we can obtain the recurrence relation for V n . Upon setting N = 2k and dividing by 
(N — 2n) we find that V n {C) is a symmetric polynomial of degree n in the variable ^ 
obeying the recurrence relation 

with V-i = 0, V\ = 1 and where 

_n(N + (-n) (N + 1 - n){n + ( - 1) 

° n ~ (2n-JV-2)' Kn = (2n - iV) ' 

As takes integer values, it is seen that when k = n, a singularity appears in the 
recurrence relation. Notwithstanding this, we proceed with the computation; the effect 
of the pole in k = n on the results is treated below. 

Introducing the monic polynomials V n (£) = , the recurrence relation becomes 



(55) 
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u. 



(N -2n)(N -2n + 2) 



The monic polynomials V n (£) can be identified with the complementary Bannai-Ito 
polynomials (CBI). 

The monic CBI polynomials [HI [16], denoted I n (x; pi, p 2 , r i> r 2), obey the recurrence 
relation 



where 



T2 



In+l(x) + (-l) n p 2 I n (x) + T n In-l(x) = xl n (x), 

n(n + pi — n + l/2)(n + pi - r 2 + l/2)(n - T\ — r 2 ) 



(2n + g){2n + g+l) 
(n + g + + pi + p 2 + + P2 ~ ri + l/2)(ra + p 2 - r 2 + 1/2) 
(2n + # + l)(2n + 5 + 2) 

and with (7 = pi + p 2 — ri — r 2 . They have the hypergeometric representation 

h n {x) = R n {x), I 2n+1 {x) = (X- p 2 )Qn(x), 



(56) 

(57a) 
(57b) 



where 



R n (x) = r] n4 F 3 

Qn{x) = L ni F 3 



—n, n + g + 1, p 2 + x, p 2 — x 
_pi + P2 + 1, P2 - ri + 1/2, p 2 - r 2 + 1/2 

—n, n + # + 2,p 2 + ;r + l,p 2 — £ + 1 
_pi + p 2 + 2, p 2 - n + 3/2, p 2 - r 2 + 3/2 



1 

;1 



(58) 
(59) 



with 



_ (pi + P2 + l)n(p2 ~ n + l/2) n (p2 ~ r 2 + 1/2) 

(n + g + l) n 
_ (pi + P2 + 2) w (p 2 - n + 3/2) n (p 2 - r 2 + 3/2) 
(n + 5 + 2) n 

and where (a) n = (a)(a + 1) • • • (a + ra — 1) is the Pochhammer symbol. 

Comparing the recurrence formulas (155]) with (I56p and (15 7p . it is seen that the 
polynomials V n (£) are monic CBI polynomials 



with 



Pn(0 = 4 (£/2;p 1 ,p 2 ,r 1 ,r 2 ) 



C - 2 2k + C 

Pi = 92 = 0, n = — - — , r 2 = 0. 



(60) 



(61) 



The parametrization (16"T|) is a special case of CBI polynomials. This case corresponds 
to the para-Krawtchouk polynomials constructed in [18J in the context of perfect state 
transfer in spin chains. 
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Since there is a singularity in the recurrence coefficients for the polynomials V n (£), 
the correspondence between the polynomials V n (£) and the CBI polynomials outlined 
above is valid only for n = 0, . . . , k and hence the Heun polynomial A(z) generates only 
the first k para-Krawtchouk polynomials. As is easily seen by induction, the recurrence 
relation (1551) generates center-symmetric polynomials V n . Hence for n > k, we have 
^n(0 = T > 2k-n{£)- Putting the preceding results together, we write 

_ (-l)M™ (k + l-n) n 

I n {U2;Pi,P2,ri,r 2 ) (62) 



for n ^ k and 



n\ (2k + (-n 
a n = a 2 k- n , n = k + 1, . . . , 2k 



for n = k + 1 • • • , 2k. 

The hypergeometric expression of the CBI polynomials (158 p provides an explicit 
formula for the coefficients a n and the coefficients b n are easily evaluated from the 
recurrence system (H5l) . Combining those results with the formulas (146p and (H71) yields 
the expansion coefficients of the eigenvectors of the operator Q in the circular basis. 
Note that these expansion coefficients only involve aj, bj with j = 0, . . . , k and hence 
only (162|) is needed. 

4-5. Eigenvectors of J 2 

To obtain the expansion coefficients of the eigenvectors of J 2 in the circular basis, it is 
necessary to relate the eigenvectors of J 2 to those of Q. This relation has been obtained 
in the previous paper [5]. In the present notation, we have 

\k, +)q = -j= (\k, +)j 2 - u k \k, -)j 2 ) , 

\ k , ~)q = -j= (\k, +)j 2 + u k \k, -)j 2 ) 
where \k, ±)j 2 are the eigenvectors of J 2 corresponding to the eigenvalues 

x ± = ^/k(k + o, 

and where the coefficient ui k is 



C-2i^k(k + Q 

Uk ~ 2kT( ' 

The inverse relation, which allows to express the eigenvectors of J 2 in terms of the known 
eigenvectors of Q reads 

1^ +)j 2 = ~j= (\k, +)q + \k, -)q) , 
\k, ~)j 2 = — 7= (| A;, +) Q - \k, -) Q ) . 
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We now consider the case fi x = \i y — //. This corresponds to a fully isotropic 2D Dunkl 
oscillator, where two "identical" parabosonic oscillators are combined. Returning to the 
system of differential equations (H9|) for the generating functions, one has 



(k + 2fi)A(z) = (k — 2/j,— zd z )B(z) + 



4/i 



X-z 1 



;B(z) 



kB(z) = (k-zd z )A(z). 

Solving for A(z), we find 

z(z 2 - X)A"(z) + (z 2 (2fi + X-2k)+2^ + 2k-X) A'(z) - Ak^zA(z) = 0. 
The solution corresponding to the initial value a = 1 is given by 



A(z) = 2 F 1 



—k, fi 
1 — k — /j, ' 



and we also have 



B(z) 



A(z) 



2kfiz 2 



k + fj, 



1 — k, 1 + fi 2 
2-k-n 



Hence in the isotropic case, the generating functions are simply the Jacobi polynomials. 
It follows from the hypergeometric generating function that 



a 2l 



i 



and 



>2n 



1 — k — /i) n rz! ' 
k-2n\ {-k) n (jj) n 



0>2n+l — 0, 



k 



[l-k- fx) n n\' 



>2n+l 



0. 



Thus it is seen that the in the fully isotropic case fi x — \i y — fi, the formulas for the 
expansion coefficients of the eigenvectors of Q simplify substantially. 



5. Diagonalization of J 2 : the iV odd case 

In this section, we obtain the expression for the eigenvectors of J 2 in the circular basis 
when N = hl + ur is an odd integer. In spirit, the computation is similar to the iV 
even case presented in the previous section. We proceed along the same lines. 



5.1. The operator Q and its simultaneous eigenvalue equation 
The operator Q is defined by 

Q = —2iJ 2 R x — /i x Ry — l^yRx — {l/2)R x R y . 
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Given the action (I24p of the reflections operators, it is seen that they have the following 
matrix representation in the circular basis B 2 '■ 



R 



1 



R x = diagOi, . . . <ri , ^ 



Using the matrix representation of J2 in the circular basis B2 given in (l3~Tj) . one finds 

A m \ 



[Gift 



with m = (JV — l)/2 and where 



/ $0 Ai A 2 

$1 Ai 



V 



A 



m— 1 



m—1 

Ai 



1/2 + i£ i(2A: + C+l)-C 

-t(2fe + c + i) -<e 1/2 -it 



A, 



-2< 


-2if 


2i£ 




2i£ 


2i£ 


-2i( 


-2t£ 



m odd 



• (63) 



/?? even 



From the block upper-triangular structure, it follows that the eigenvalues of Q are 



u+ = 2k + C + 3/2, 



/A. 



-(2A + C + 1/2), 



for = 0, . . . , m. Although a different labeling has been used, it is directly checked 
that the eigenvalues of Q are the same for the N even and iV odd case, except for the 
additional one. We denote the eigenvectors of Q corresponding to the eigenvalues 
by \k, ±)g and define their expansion in the circular basis by 

m m 



<T=± 



for k — 0, . . . , m and where the vectors \£, ±) are the vectors of the circular basis $2. 

We shall once again study the simultaneous eigenvalue equation for the operator 
Q. We define the matrix of eigenvectors 



W 



fV 00 V01 
V n 



Vim 
VmmJ 



where 



vj{k) uj(k) J ' 
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The simultaneous eigenvalue equation for the matrix [Q]g 2 reads 

W-L=[Q] B2 W, (64) 

with 

L = diag(A , • • • , A m ), A fc = ^ . (65) 

As in section 4, the simultaneous equation (164)) will be shown to be equivalent to a 
system of recurrence relation for the components Ue(k) ± , v l {k) of the eigenvectors of Q 
in the circular basis. 



5.2. Recurrence relations 

We now construct the recurrence systems for the component of the eigenvectors of the 
operator Q. For k = 0, . . . ,m and £ = 0, . . . ,m, the simultaneous equation (164|) takes 
the form 

k-i 

V ek A k = $ e V ek + ^J v jk- (66) 
i=i 



5.2.1. The v k eigenvalue sector 

Let us begin by considering the eigenvectors \k, — )q of Q corresponding to the eigenvalue 
vj~ and their expansion in the circular basis 

k 

i*.->e=X>f(*)iM- 

e=o 

It is directly seen from (1661) . (163|) and (165)) that we have 

k k 

[2k + 1 + C + iftvf = [f - *'(2* + C + - 2^ ^ (-l) j ~% - 2*^, Bj, (67a) 

j=M-l j=£+l 
fc k 

[2k + l + (- i£]vj = {£ + i(2£ + C + 1)K+ + 2i/x a ^ (-l) j ~ e Aj - 2ifi y ^ S,-, (67b) 

j=t+i j=e+i 

where we have defined Aj = + v J and Bj = vj — vj and where the explicit 
dependence on k of the components vfijz) has been dropped for notational convenience. 
The recurrence system (1671) is "reversed". The terminating conditions are at £ = k. In 
this case, ( ISTj) becomes 

[2k + 1 + C + %W = K - *(2* + C + 1)K, 
[2fc + 1 + C - = [£ + i(2As + C + 1)K + - 
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Choosing v7 = 1, we obtain the terminating conditions 

v k = ~h % = 1. 

Upon taking 

A t = a A £} B e = O B £ , 

where ckq = (1 — i) and /3q = (1 + i), it is easily seen that Ai and are real and that 
the system ( 1671) is equivalent to 

k 

[k + fiy + l/2]A e = (£ + fi y + l/2)B e + 2u y % 

j=e+i 

k 

[k + fi x + l/2]B e = {£ + fi x + 1/2) A t + 2fi x 

j=l+l 

with the terminating conditions A^ = 1 and B^ = 1. Introducing the reversed 
components a n = Ak- n , b n = Bk- n , we obtain the system 

n-l 

[k + Hy + l/2]a„ = [fc - n + ^ + l/2]6 n + 2/^ (68a) 

n-l 

[A; + + l/2]b n = [k-n + fi x + l/2]a n + 2fi x {-l) n ^(-1)^, (68b) 

i=o 

with the initial conditions So = 1 and bo — 1. Taking into account all the preceding 
transformations, we have 

v e \ K ) — 2 ' v i ^ ' ~ 2 ' 

where a^, bi are the unique solutions to the recurrence system (168]) . 

5.2.2. The u£ eigenvalue sector 

Let us now consider the eigenvectors \k, +)q corresponding to the eigenvalue v£ with 
expansion 

k 

|*,+>Q = X>J(*)I4">. 
e=o 

in the circular basis. Proceeding along the same lines as in the previous computation, 
we find that the terminating condition are of the form 

+ _ i(2fc + l + C + »Q __, 
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and that the components are given by 



u e (k) = , uj{k) = , (70) 



where 



_ (l + i)(2fc + l + 2/j y ) _ (1 - i)(2k + 1 + 2fi x ) 
70 ~ 2k + l + £-i£, ' e °~ 2k + l + (-i£, ' 

5.3. Generating functions and Heun polynomials 

As is seen from the formulas (169 p and (|70|) . the main part of the components of the 
eigenvectors of Q in the circular basis is given by the solutions to the recurrence system 
(j68]). As in section 4 for the N even case, we bring the ordinary generating functions 

A(2)=X)V, B(z)=Y}nZ n . 
n n 

Upon using the identities (I48p . we obtain from (168]) the associated differential system 
(A; + fi y + 1/2)A(2) = (k — fly + 1/2 - zd z )B(z) + -£!-B(z), (71a) 
(fc + ^ + 1/2)50*) = (k - fi x + 1/2 - ^)l(z) + P^-A(z). (71b) 

1 + -2 

By direct substitution, we find that the generating functions are expressed in terms of 
the Heun functions 

A(z) = (1 + z)H£(a, q A ; oi A , (3a, 1a, 6a, z), (72a) 
B(z) = (1 - z)H£(a, q B ; a B , (3 B , 1b, 6 b , -z), (72b) 

where 

a = -1, q A = 2k(p y -fi x - 1), ol a = -2k, (73a) 
/3a = + + 1, 7a = ~2k - fu x - fu y , 8 A = 2fi y . (73b) 

and where the parameters qB,'~'8 B are obtained from ( 173)) by the transformation 
/■tec ^ AV The form of the parameters involved in the Heun functions show that once 
again one has a truncation at degree 2k + 1 and hence the Heun functions appearing 
in (172"]) are in fact Heun polynomials. The generating functions A (z) and B(z) are 
polynomials of degree 2k + 1. 



5.^. Expansion of Heun polynomials in complementary Bannai-Ito polynomials 

The expansion of the Heun polynomials can be obtained using the associated three-term 
recurrence relation. Since the expansion coefficients of B(z) and A(z) are related by the 
simple relation z -H- —z and fi x -H- //j,, we shall focus on the expansion of A(z). 
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We examine the expansion of the Heun polynomial appearing in (I72ap 

H£(a, q A , a A , (3 A , ^ A , 5 a , z) = ^ ?V& r , (74) 

r 

where the parameters are given in (173]) . Using the recurrence coefficients given in (154|) . 
one finds that the expansion coefficients V n are polynomials of degree r in £ = fi x — fi y 
that obey the recurrence relation 

*r+lPr+l(0 + «r^r-l(0 = (1+0^(0, (75) 

with V-i — 0, Vq — 1 and where 

(r + l)(iV + C-r) (N + l-r)(r + () 

ar+1 ~ 2r~^N ' Kr ~ 2r~^N ' 

where we have taken N = 2k. Introducing the monic polynomials V r (£) = J^..^ , we 
obtain 

Vr+l(0 + VrVr-1^) = (1 + 0^(0 (76) 

where 

r(N + 1 - r)(N + ( - r - l)(r + () 

11 — 

(N-2r)(N-2r + 2) 

Upon comparing (176]) with the recurrence coefficients for the CBI polynomials given in 
(I57jl . it is directly seen that the polynomials defined by the recurrence (1761) correspond 
to CBI polynomials with the parametrization 

C-l 2k + Q + l (77 , 

Hence, when r ^ k, we have 

?r(fl = ( ~ 1 J 4 \ 2 f + + / + 7-r)/ - ((1 + 0/2;Pl ' P2 ' ri ' r2) ' r °' (78) 

where pi, p 2 , ?"i and r2 are given by ( 1771) and J n (x; pi, p2, r±, r 2 ) are the complementary 
Bannai-Ito polynomials. As can be seen by the recurrence relation ( 175]) . the expansion 
coefficients of the Heun function (1741 truncate at order 2k + 1 and are center-symmetric. 
Hence, for r > k, we have 

Pr(0 = ?Wl-r(0: ^ > k. (79) 

Taking into account the relation between A(z) and B(z), the expansion coefficients of 
the Heun function appearing in (]72b|) . denoted by T n {C)i are easily seen to be 

^(0 = ^ (2 f + + c 1 + ~ r _j r 1) / r((l-0/2;Pi,P2,r 1 ,r 2 ), r^2k, (80) 
T r (0 = W-r(0, > 2*. (81) 



The Dunkl oscillator in the plane II 27 
Collecting all the previous results, we write 

a n = V n (Z)+Vn-i(0, (82a) 

&n = r„(o + r„-i(o, (82b) 

where V-i = 0, 71i = and V n (Q, %(£) are given by ([78]), ([79]), (EDJ and (|gl)). 
5. 5. Eigenvectors of J 2 

To obtain the expansion of the eigenvectors of J 2 in the circular basis, one must relate 
the eigenvectors of the operator Q to the eigenvectors of J 2 . This relation has been 
obtained in the previous paper |5J. We have 

\k, ±)q = -j* (\k, +) j2 =F v k \k, -) Ja ) 

where \k, ±)j 2 are the eigenvectors of J 2 corresponding to the eigenvalues 



Xf = ±^{k + fx x + l/2){k + fi y + 1/2), k = 0, . . . , m, 

and where 

+ 2V(fe + ^ + V2)(fe + % + 1/2) 



The inverse relation reads 



2fc+C+l 



\k,+)j a = ^={\K+)a+\k,-) Q ), (83a) 

\k, -). h = -^L flfc, +) Q - -) c ) (83b) 



V2 



Using the relations f l83l) . the results (l8~2~l) . f l70l) and f l69l) . one has an explicit expression 
for the expansion of the eigenvectors of J 2 in the circular basis for the case iV odd. 

5.6. The fully isotropic case : —1 J acobi polynomials 

We consider again the case /j, x — \i y — \i which corresponds to the fully isotropic 
Dunkl oscillator, where two independent identical parabosonic oscillators are combined. 
We return to the system of differential equations for the generating functions of the 
components of the eigenvectors of Q given in ( 17TT) . When fi x = \i y = fi, one has 

(k + fji+ 1/2)2(2) = (k - fi + 1/2 - zd z )B{z) + -^-B(z), 

1 — z 

{k + fi + l/2)B(z) = (k - fi + 1/2 - zd z )A(z) + -^—A(z). 
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It is easily seen from the above formulas that A(z) = B(—z). Hence the generating 
function A(z) satisfies the differential equation 



(k + fj,+ 1/2) A{z) = (Jfe - ii + 1/2- zd z )A{-z) + 



2/i 



1-z 



A(-z), 



which may be cast in the form of an eigenvalue equation 

LA(z) = 4:/j,A(z), 

where 

2£; + 1 + 2/i 



L = 2(1 — z)d z R + 



;-2fc-l + 2 / u) 



(I-R), 



where Rf(z) = /(— z). It is recognized that the operator L is a special case of the 
defining operator of the little —1 Jacobi polynomials [17J. 

The little —1 Jacobi polynomials, denoted by P~ l (x), obey the eigenvalue equation 



QP n ( X ) = \ n Pn{x) 



A, 



— 2n n even 

2(a + (3 + n + l) n odd 



where 



ft = 2(1 - x)d x R + (a + p + 1 - a/x)(I - R). 



Comparing the operators Q and L, it is seen that the generating function A(z) 
corresponds to a —1 Jacobi polynomial of degree n = 2k + 1 with parameters 



a = —2k — 2/^ — 1, 



/3 = 4/i- 1. 



Using this identification and the explicit formula for the little —1 Jacobi polynomials 
derived in |17j . we obtain 



A{z) = 2 F 1 



-fi — fc' 



+ ~ — : — 2-ri 



-k,/J, + l 

-k — n + 1 ' 



k + /i 

This directly yields the following result for the recurrence coefficients a 

'-k) n (n) n „ fl (-k) n (/J, + l) n 



n\{—ji — k) n ' ' k + iin\(—[i — k + I) 
Using the symmetry B(z) = A{— z), we also obtain 



J2n 



1 



n\(-fi - k) n 



-fi (-k) n (fi + l) n 

k + jj, n\(—/i — k + l) n 



Once again, the components of the eigenvectors of Q drastically simplify in the isotropic 
case /j^ x = /Xy = /x. Moreover, the preceding computations entail a relation between a 
special case of Heun polynomials and the little —1 Jacobi polynomials. 
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6. Representations of sd{2) in the J 2 eigenbasis 

We now investigate the representation space in which the operator J2 is diagonal. The 
matrix elements of the generators of the Schwinger-Dunkl algebra will be derived using 
the defining relations of sd(2), which read 

{J 1 ,R Xt } = 0, {J 2 ,R Xl } = 0, [J 3 ,R Xi } = 0, 

[J2,J3]=iJl, [J3,Jl]=iJ2, 
[Jl, J 2 ] = i[J 3 + J^xRx + VyRy) - U{ii x R x - II y Ry)/2], 

where R? x = R? y = I. It will prove convenient to treat the even and odd dimensional 
representations separately. 

6.1. The N odd case 

We first consider the case where N is odd. The representation space C is spanned in 
this case by the basis vectors \k, ±) with k G {0, . . . , m} on which the generator J 2 acts 
in a diagonal fashion 

J 2 \k,±) = X^\k,±), k = 0,...,m, 

where m = (N — l)/2 and where the eigenvalues of J 2 , derived in section 3, are given 
by 

A± = ±^J(k + /2 x + l/2)(fc + /Xj, + l/2), k = 0,...,m. 

Since R x , R y anti-commute with J 2 and given that R x R y is central in the algebra sd(2), 
we can take 

R x \k,±) = e\k,T), R y \k, ±) = \k, =f), 

where e = ±1. Here we choose e = —1, which corresponds to the representation 
encountered in the model. In the basis {|0, +), |0, — ), • • • , \m, +), \m, — )}, the matrices 
representing the involutions R x , R y have the form 

R y = -R x = diagOi, . . . , at), o x = 

which is identical to their action in the circular basis. The Hamiltonian "H has the action 

H\k,±) = (N + n x + n y + l)\k,±). 

The action of the operator J3 on this representation space can be derived by imposing 
the commutation relation using ([6]) to define J\. The action of J 3 on the basis \k, ±) 
is taken to be 

J 3 \k,+)=J2M%\j,±), Mk,-) = J2Nf k \j,±). 

3=0 j 
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With these definitions, it is easily seen that the commutation relation ([2]) is equivalent 
to the following system of relations 

C(iV + 1 + C) 



2(A+A+) + (A+) 2 



[^) 2 -2(A-A+) + (AT)2 



P fc ) 2 



2(A+A^) + (A+)< 



[(Ar) 2 -2(A-Ar) + (AT)^ 



jk 

- mi 



(Pv-Hx)N7 k + 8 jh : 9 

C(iv + i + C) 



Hx)M. 



(84a) 
(84b) 

(84c) 
(84d) 



where the relations (I84a|) . (I84b|) were obtained by acting on +) and the relations 
(I84cp . p4d|) by acting on \k, — }. It follows directly from the solution of the system (jS] 
that J3 acts in a six-diagonal fashion on the eigenbasis of J 2 . For j = k, we obtain 

tt(N + ( + l) . C (at + C + i) 



A'/,- 



2(2fe + C)(2/e + C + 2) 

C(iv + C + i) 



A', 



A'/,- 



fcfe 



A r , 



2(2fc + C)(2/c + C + 2) ' 

ec(iv + c + i) 



kk 



(85a) 
(85b) 



2(2fc + C)(2/c + C + 2)' 
or j = k — I with I > 1, only the trivial 



2(2fc + C)(2/c + C + 2 

where £ = /i x - /j, y , £ = to + to- For 3 = k + 
solution occurs, so the matrix representing J3 in the eigenbasis of J 2 is block tridiagonal 
with all block 2x2. Hence it acts in a six-diagonal fashion on the eigenbasis of J 2 . Using 
the commutation relations and the system (I84p . it is possible to obtain an expression 
for the matrix elements of J3 which involves a set of arbitrary non-zero parameters {f3 n } 
for n — 0, . . . , m. After considerable algebra, one finds that the matrix J 3 has the form 



/Co 
Do 



■h 



where the blocks are given by 



U-i 



Co 



D 



771—1 



u k = p k 



M 



k-lk 
1 



M, 



k-lk 



Ck 



kk 



M kk N kk 



\ 



Ufn 
C m ) 



D k = B 



k+l 



( M k + lk 

W+ik 



The matrix elements of the central blocks are given by ( 1851) . The components of the 
upper blocks U k have the form 

' fl - 4(* + u x )(k + to)~4{(k + u x - l/2) 2 {k + u y - 1/2) 2 } 1/2 1 



M 



k-lk 



2£ 



The matrix elements of the lower blocks have the form 

£(k + l){2k-N + l)(k + l + ()(2k + 2( + N + 3) 



K+ik 



4(2fc + C + 2)(2A: + C + l) 



k+lk 



E k M, 



k+lk- 



where 



Eu — — 
2£ 



1 - 4(fc + u x + l)(fc + ^ y + 1) + A{{k + u x + l/2) 2 (k + to + 1/2)2} 



1/2 



We note that these matrix elements are valid for \i x 7^ \i y - m the latter case, the form 
of the spectrum of J 2 changes and the computation has to be redone from the start. 
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We now consider the N even case. The representation space C is spanned in this case 
by the basis vectors |0, — ) and \k, ±) with k — 1, . . . , m on which the operator J 2 acts 
in a diagonal fashion 

J 2 \k,±) = \ k \k,±), k = 0,...,m, 

where the eigenvalues of J 2 , determined in section 3, are given by the formula 

\± = ±y/k(k + (), fc = 0,...,m. 

Note that the eigenvalue A is non-degenerate. In this representation, we choose the 
following action for the involutions R x , R y : 

R Xi \0,-) = \0,-), R Xi \k,±) = \h,^), 
and hence the reflections have the matrix representation 



i? a . = i2 w = diag(l,<Ti,...,<Ti), o-i 



1 

1 



The central element (Hamiltonian) H has the familiar action 

n\k, ±) = (N + fji x + Hy + 1) | At, ±). 

Following the same steps as in (6.1), a direct computation shows that the in this case 
J3 has the matrix representation 



■h 



( CO Ul 

do Ci U 2 
Di C 2 

V 



\ 



'•■ u m 

Dm-l Cm/ 



The special 2x1 and lxl blocks are given by 

co = — 2 ^ - ^ — , ui = [ai cti), d = a x [w N w N ) , 

with 

(iV/2)(l + 2u x )(l + 2n y ){N/2 + C + 1) 



The 2x2 blocks have the form 



2(1 + C) 2 (2 + C) 



m-aJ 1 N t~A C h -( M ^ N ^A D h -a~^ ( M ^ 
Uk ~ ak W_ lk 1 J' Ck ~W k M+J' Dk ~ a ^W +lk M+ + J- 
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where 

+ _ ttiN + C + l) + _ -fliV + C + 1) 



kk 2(2fc- l + 0(2fc + l + C)' fefe 2(2fc- 1 + C)(2fc + 1 + C)' 
AT+_ lfe = C 1 {C + 2(fc - l)(fc + C) ~ 2[(fc - l) 2 (fc - 1 + C) 2 ] 1/2 } , 

+ _ (N/2 - k){N/2 + fc + 1 + Q(2fc + 1 + 2fi x )(2k + 1 + 2/x„) {C + 2fc(fc + ( + !) + 2[(fc) 2 (fc + Q 2 ] 1/2 } 

4(2fc + l + C)(2fc + C) 3 
+ _ C(N/2 - k)(N/2 + k + C + l)(2fc + 1 + 2^)(2fc + 1 + 2/V) 



4(2fc + C + l)(2fc + C) 3 

The parameters of the sequence are arbitrary but non-zero; they could be fixed, for 
example, by examining the action of J2 on the eigenstates of the 2D Dunkl oscillator in 
the polar coordinate representation. We have thus obtained the action of the operator 
J3 on the eigenstates of J 2- Recall that J3 is the symmetry operator associated to the 
separation of variables in Cartesian coordinates and J 2 is the symmetry associated to 
the separation of variables in polar coordinates. 



7. Conclusion 



We have investigated the finite-dimensional irreducible representations of the Schwinger- 
Dunkl algebra sd(2), which is the symmetry algebra of the two-dimensional Dunkl 
oscillator in the plane. The action of the symmetry generators in the representations 
were obtained in three different bases. In the Cartesian basis, the symmetry generator 
J 3 associated to separation of variables in Cartesian coordinates is diagonal, and the 
symmetry J2 is tridiagonal. In the circular basis, the operator J3 acts in a three- 
diagonal fashion and J2 has a block upper-triangular structure with all blocks 2x2. The 
eigenvalues of J2 can be evaluated algebraically in the circular basis and the expansion 
coefficients for the eigenvectors of J2 in this basis are generated by Heun polynomials 
and are expressed in terms of the para-Krawtchouk polynomials. Finally, it was shown 
that in the eigenbasis of J2, the operator J3 acts in a block tridiagonal fashion with all 
blocks 2x2, that is, that J 3 is six-diagonal. 

It has been seen that the Dunkl oscillator model is superintegrable and closely 
related to the —1 orthogonal polynomials of the Bannai-Ito scheme. In this connection, 
the study of the 3D Dunkl oscillator model and the singular 2D Dunkl oscillator 
could also provide additional insight in the physical interpretation of the orthogonal 
polynomials of the Bannai-Ito scheme. 
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